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Abstract. We study the evolution of wormhole geometries under Ricci flow 
using numerical methods. Depending on values of initial data parameters, 
wormhole throats either pinch off or evolve to a monotonically growing state. 
The transition between these two behaviors exhibits a from of critical phenomena 
reminiscent of that observed in gravitational collapse. Similar results are obtained 
for initial data that describe space bubbles attached to asymptotically fiat regions. 
Our numerical methods are applicable to "matter-coupled" Ricci flows derived 
from conformal invariance in string theory. 



1. Introduction 

An analog of the diffusion equation for geometries is the so-called Ricci flow. This is a 
type of parabolic partial differential equation that homogenizes geometry in a manner 
similar to the homogenization of density or heat produced by the diffusion equation. 

The Ricci flow equation was first studied in the early 1980s. It arose nearly 
simultaneously in mathematics [1], in the study of the classification of 3-geometries, 
and in string theory physics [2j [3] from the requirement that a quantized string retain 
conformal invariance. Certain modifications of the purely geometric flow equation by 
"matter terms" motivated by string theory were central to Perelman's [4] work on the 
Poincare conjecture. 

A goal of analytic work on the flow equations is to study the approach to 
singularities, and to understand the so-called "ancient solutions," which are the non- 
linear analogs of the heat kernel of the diffusion equation. Much of the work in 
this area has focused on the flow of compact geometries [S], with relatively little on 
non-compact cases such as the asymptotically flat or anti-deSitter geometries. In the 
asymptotically flat case for dimension n > 3, it has been demonstrated recently that 
the flow exists and remains asymptotically flat for an interval of time, with the mass 
remaining constant [6l[7]. For the rotationally symmetric case, Ref. [6] also shows 
that the flow is eternal for data containing no minimal surfaces. 

The purpose of this paper is to investigate numerically some examples of the 
flow of asymptotically flat geometries, and to describe techniques that are easy to 
generalize to cases where matter flow equations are included. 

Previous numerical work on Ricci flows focused on the (compact) case of a 3- 
sphere that is corseted (i.e., a surface formed by uniformly shrinking the equator of a 
3-sphere), where it was found that the manifold flows to a round 3-sphere or a neck 
pinch depending on the degree of corseting i8j . The flow of horizon area and Hawking 
mass has also been studied |9] for initial data that is viewed as time symmetric data 



Ricci flows, wormholes and critical phenomena 



2 



for the Einstein constraint equations. Numerical work pertaining to string theory aims 
at obtaining new black hole solutions as fixed points [TD] . 

The focus of the present work is the geometric flow of wormhole geometries. 
We study several classes of initial data and find two possible outcomes: the wormhole 
"pinches off" leaving two disjoint asymptotically flat regions, or it expands indefinitely. 
We also find that there is a type of critical behaviour that separates these extreme 
cases. The numerical methods we utilize are stable for long time flows, and have 
potential applicability to the largely unexplored geometry-matter cases. 

2. Flow equations 

The Hamilton-DeTurck Ricci flow is defined by the equation [TJ [TT| 

dt9ab^-2Rab + 2V^aVb). (f) 

For our work, we take gab = gab{t,x'^) to be the 3-metric on a manifold St with 
Euclidean signature and Rab is the associated Ricci tensor. The parameter t labels 
different 3-geometries along the flow and should not be confused with the intrinsic 
3-dimensional coordinates x'^. Finally, V"" is a DeTurck vector field that generates 
diffeomorphisms along the flow. The DeTurck field essentially represents the freedom 
to change 3-dimensional coordinates as the flow progresses. It is useful to note that t 
has dimensions of [length]'^ and has dimensions of [length]"^. 

We assume the 3-manifolds are spherically symmetric, for which a general 
metric ansatz is 

ds^ = e^^^'^P^dp" + R\t,p)d^f]. (2) 

Spherical symmetry implies the following form for the DeTurck field: 

V'^ = V{t,p)dp. (3) 

When ^ and ([3]) are put into (H]) we obtain two independent flow equations for 
dtX{t,p) and dtR{t,p). There is no dynamical equation for V{t,p), so the system 
is underdetermined, which reflects the coordinate freedom embodied by the DeTurck 
field. To close the system, we need to fix coordinate gauges. There are three choices 
we have investigated: 

(i) V{t,p) = 0; 

(ii) R{t,p) = p (conformally flat gauge); or 

(iii) X{t,p) = (areal radius gauge). 

Enforcement of condition ([u]) or (pS)) in the flow equations ([T]) results in a single 
dynamical equation for X{t,p) or R(t,p), respectively, and an equation of constraint 
for V{t,p). The choice Q gives two dynamical equations for both metric functions. 

For our purposes, a wormhole geometry is defined by a metric ([2|) that has two 
asymptotically flat regions as p ^ ±c». The combination e-^'^'*''''^ R{t, p) should be 
non-zero for all p, and the wormhole throat (or throats) are located at local minima 
p = Pth of this functioijl 

dp{e^R)\ ^0, dl{e^R)\ > 0. (4) 
'\p=Pth ' ''p=pth ^ ' 

% If we interpret the 3-geometries St as being embedded in 4-dimensional Lorentzian spacetimes as 
moments of time symmetry, each wormhole throat represents an apparent horizon. However, if the 
Ef are embedded with non-zero extrinsic curvature, such an interpretation is not possible. 
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The definition of the wormhole throat suggests a convenient coordinate gauge choice 
is given by case (jm]) above, with hne element 

ds^ = dp^ + R^{t,p)dn^. (5) 

In this gauge the area of 2-spheres of constant p is 47ri?^(<, p), hence the name "areal 
radius" gauge. The condition X = reduces the flow equations ([T]) to 

jdpR)^ 1 
R R 
dlR 

dpV^^2^. (66) 

As mentioned above, the latter is a constraint equation. These equations must be 
supplemented by boundary conditions and initial data. In the areal radius gauge, 
asymptotic flatness requires that 

R^\p\+0{h\n{\p\/h)). (7) 



dtR = dlR + '-^^ -- + VdpR, (6a) 



where li and I2 are two length scales. Putting this asymptotic limit into (j6a|l and (|6bj) . 
we obtain 

V^O{\p\-'). (8) 

The asymptotic behaviour of R suggests that we search for solutions that are even 
functions of p. We then only need to consider the interval p G [0, 00), and we have the 
additional boundary condition 

The boundary conditions imply that there is always a local extrema of R at p — 0. If 
this is a minima, then there is a throat at p = 0. To complete the specification of the 
flow, we need to give initial data for R at t = 0. The various classes of initial data 
considered in this work are described in 21 ^ 

Finally, we note that it is possible to embed the 3-geometries ([5]) in 4-dimcnsional 
flat space 

ds^ = dR^ + dZ^ + R^dn^, (10) 
via the parametric equations 



R^Rit,p), Z = Zit,p), dpZ=^l-{dpRY. (11) 

These embeddings are only possible for areal radius functions satisfying \dpR\ < 1. 
We use these as a tool for the visualization of Ricci flow initial data in Figs. [T]and[H 



3. Numerical method 

We solve the equations (|6ap and (j6b|) numerically using finite difference methods. 
The computational domain is defined by p € [0,/5niax] and t € [0,tniax]- The spatial 
and temporal intervals are discretized into elements of size Sp and 6t, respectively. 
The notation Aj refers to the value of the quantity A{t, p) at the i^^ spatial and j^^ 
temporal node. The Neumann boundary condition ([9]) is imposed explicitly, while the 
asymptotic boundary conditions ([7]) and ^ are replaced with the Dirichlet conditions 

R{t, Pmax) = i?max, V{t, p^ax) = 0. (12) 
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Simple Euler (SE) 



Modified Dufort-Frankel (MDF) 



dpA 
dpR in eq. 
dpR in eq. ^ 

dtR 
accuracy 



{Ai^,-Ai_,)/{2Sp) 



■i+l 



{Rr~Ri)/St 
0{5t,5p^) 



2Rl 
2R^ 



^i-i)/(2M 

-1 



{Rl^^-R^^-Rf^+R^_^/5p^ 
m^,-2R^+^_^/5(P 

{Ri+'~Rr')/i2St) 



Table 1. Summary of the simple Euler (SE) and modified Dufort-Frankel (MDF) 
finite difference derivative approximations. In the top line, A stands for either R 
or V. Notice that for the MDF method, the second-order spatial derivative of R 
is handled differently in the flow 116 all and constraint 1 I6&I I equations. 



Rmux is fixed by the asymptotically flat initial data. With this, the former condition 
efTectively gives dR/dp ^ 1 for large p in our simulations, while the latter leads to 
a non-zero value of V at the origin p = ol^ For pmax sufficiently large compared to 
imax, the simulation results are insensitive to the position of the outer boundary. In 
practice, our simulations are performed with i?max ~ Pmax > 200. 

The two types of finite differencing schemes used to approximate derivatives 
are summarized in Table [T] When applied to linear parabolic systems, the simple 
Euler (SE) method is known to be conditionally stable when 5t/5p^ < 1. Conversely, 
the modified Dufort-Frankel (MDF) method is known to be unconditionally stable. 
Via direct experimentation, we have found that these conclusions also hold for the 
nonlinear problem given by equations (j6ap and (|6&P , which would naively suggest that 
the MDF method is vastly superior to the SE method due to its excellent stability. 
However, there is a non-trivial price to be paid for this attractive feature. One can show 
that the numerical results obtained with the MDF approximation actually converge 
to solutions of 

m = dlR + _ ^ + vd,R (I) ' dfR, (13) 

rather than the solutions of the original flow equation (|6ap fj| In other words, the MDF 
method solves a modified version of the original system of PDEs. In order to minimize 
the discrepancy between MDF results and the "true" solutions of (|6aP , one must make 
the ratio e = St/Sp as small as is computationally feasible. 

The numerical results presented in this paper are obtained using the MDF method 
with e = 10"''. For several individual cases, we have checked that the simulation results 
are insensitive to the particular choice of 6t and 6p provided that e < lO^'^, and that 
the MDF and SE schemes give virtually identical answers for e = 10^"*. 

§ This means that data symmetric about the origin evolves asymmetrically. Since our computational 
domain is p £ [0,pmax], the other side of the wormhole (ie. p G [— Pmax, 0]) can be glued continuously 
using the junction conditions with a jump in V at p = 0, which represents the freedom to choose 
different coordinates on each side of the join. 

II The origin of the anomalous time derivative on the r.h.s. of (I13I I is the MDF approximation of d^R 

in I l6al l. Essentially, the MDF prescription approximates Rj by its temporal average (R^^^ +-R^~^)/2. 
The use of nonlocal time data to resolve the spatial derivative gives rise to the extra dfR term. 
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a = 1 (Einstein-Rosen bridge) a = 2 (Ellis wormliole) 



Figure 1. Embedding diagrams of the Einstein- Rosen and Ellis wormholes in 
fiat space (we have set 6 = it/2 to obtain a 3-dimensional plot). 

4. Evolution of Morris-Thorne wormhole geometries 

The Morris-Thorne wormholes 12] have spatial sections 

ds^ = -i-dr2 + r^dn^, f{r) - 1 - (14) 
/(r) r 

In this section, we will use these metrics as initial data for the Ricci flow. At the 
initial time t = 0, this 3-metric can be transformed into the areal radius gauge ([S]) by 
the coordinate transformation r = i?(0, p) with 

9pi?(0,p)- VTMM), (15) 

where i?(0, p) is the initial profile for the areal radius metric function. We will 
concentrate on a subset of this class defined by 

b(r) = ^, m = (16) 

where rg is a length scale and a is a dimensionless parameter (a similar family of 
wormholes was defined in 13J). The spatial sections of the Schwarzschild metric 
(known as the Einstein- Rosen bridge) correspond to a = 1. Another special case 
is the Ellis wormhole with a = 2 [14], which has R{0,p) = vT^+Tq. Embedding 
diagrams of both types of wormholes appear in Figure [1] For Ricci flow of these 
wormholes, it is useful to work with dimensionless quantities defined by 

p^p/ro, i^t/rl R^R/ro, V ^ roV. (17) 

When the hatted-variables are substituted into the fiow equations the tq length scale 
drops out of the problem; i.e., our simulation results are essentially independent of tq 
up to trivial scalings. 

In Figure [21 we show the results of our simulations for the evolution of the areal 
radius function R{t, p) for the Einstein-Rosen and Ellis cases. One can see in Figure 
[1] that the initial data for the two scenarios are visually quite similar. However, the 
evolution of R is entirely different: in the Ellis case the areal radius tends to diverge 
for all /9 as i — > CX3, while for the Einstein- Rosen bridge R is everywhere decreasing 
and the wormhole throat pinches off in finite time. Due to computing limitations, we 
cannot determine the ultimate fate of the Ellis initial data, but it appears to be a 
warped cylinder with a monotonically increasing radius. 

In Figure[31 we plot the evolution of the throat radius Rth{t) = R{t, 0) for various 
values of a. We find that the wormhole throat pinches off for initial data with a 




Figure 2. Evolution of the areal radius function R{t, p) under tlie Ricci flow for 
Morris-Thorne wormhole initial data corresponding to a = 1 and a = 2. Note 
that our boundary conditions ensure that the wormhole throat is always at p = 
for these simulations. 




flow parameter I/tq flow parameter log£/ry 



Figure 3. Ricci flow evolution of the throat radius -Rth(t) = R{t,0) for various 
Morris-Thorne wormhole initial data. 



less than a critical value ctcrit- For cases with a > ctcritj the throat approaches 
an ever-expanding state after a brief initial period of contraction. Numerically, we 
have determined the value of the critical parameter to be approximately 1.259. For 
a > Qicrit, the throat radius appears to grow like ^/t at late times. 

5. Evolution of "bubble" geometries 

A geometrically interesting question is: What does the Ricci flow evolution of a bubble 
attached to asymptotically flat space look like? In the areal radius gauge, initial data 
for this kind of problem would roughly be 

Tjfn \ Jsmp/ro for < p < po 
[ p/ro for p> po 

where po < irro is some transition radius. Notice that this initial data has R = and 
dpR ^ at /9 = 0, so it violates the boundary conditions that our code was designed 
to deal with. We therefore consider the related problem of a bubble attached to two 




Figure 4. Embedding diagrams of 3-geometries representing a bubble connecting 
two asymptotically flat regions 



asymptotically flat regions, and leave initial data of the form (fTSl) for future work. We 
take initial data for the current problem to be 

(ro{l-pyrf,) for|p|<po 

R{0,p) = { I ■ (19) 

[\J{\P\- Pi? + pI for|p|>po 

Here, pi and p2 are selected such that both R and dpR are continuous at the transition 
points p = ±po- Defining the dimensionless parameter j3 = po/ro, we show embedding 
diagrams for two examples of this type of initial data in Figure ID Note that there are 
two wormhole throats (local minima of R) in this case at p = ±pi, while there is a 
local maxima of i? at p = 0. This maxima corresponds to the equator of our bubble, 
and it is easy to show that the ratio of the throat radii to the equator radius is 



i?th/^cq = P2/ro - (1 - /3')v/l - 4/32. (20) 

Hence, in order to have sensible initial data with i? > we are obliged to take 
0<P< 1/2. 

Making use of the dimensionless quantities ([T7)) , we have simulated the Ricci flow 
evolution of various types of bubble geometries. As in 21 above, we take R to be an 
even function of p and impose the boundary condition that dpR = for p = 0. We find 
that for Rth/Req initially small, the bubbles disappear and the geometry "pinches-off" 
into a pair of disjoint asymptotically flat regions. For i?th/-Req initially close to unity, 
the bubbles also disappear and the areal radius blows up, somewhat similar to the 
ct > Qfcrit cases in gl Explicit examples of these behaviours are shown in Figure O 



6. Conclusions 



We have observed that the Ricci flow of certain families of rotationally symmetric 
wormhole geometries exhibit a form of critical behaviour: wormhole throats pinch off 
or expand forever depending on the values of initial data parameters. The expanding 
solutions are perhaps surprising at first sight, but are nevertheless understandable in 
that they appear to converge to an infinite volume cylindrical attractor. 

This work complements both the analytic work on the flow of asymptotically 
flat geometries [6] and the numerical work on corsetted 3-spheres [8]. It would be 
of interest to see if the neck pinching seen in our simulations of wormhole throats 




Figure 5. The evolution of two bubble geometries under the Ricci flow. 



may be modeled by the Bryant steady solitons, as has been observed for rotationaUy 
symmetric geometries on [15] . 

The methods we have used have appUcabiUty to a variety of situations for 
relevance to string theory, where the renormalization group flows give matter coupled 
modifications of Ricci flow. The interest here is in discovering new fixed points of 
physical relevance. For example, in special cases with compact target spaces and a 
tachyon field, there is evidence that properties of the tachyon potential are linked to 
the existence of non-trivial fixed points f6j. ft is of interest to extend such cases to 
non-compact geometries with gauge or 2-form fields, and consider questions such as 
the possibility of stabilization of wormholes geometries. 
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